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ABSTRACT 



Reoorted are the findings of the Uspekhi experiment 
in creating a labeled machine file, as well as sample products of 
this system - an article from a scientific journal and an index page. 
Production cost tables are presented for the machine file, primary 
journals, and journal indexes. Comparisons were made between the 1965 
predicted costs and the actual 1969 costs. Despite the greater 1969 
production cost, the authors anticipated cost reductions with th.e 
implementation of advanced equipment; such as Cathode-ray tube (CRT) 
equipment, on-line PDP-8 editing equipment, adequate Photon disc and 
translation bridge from Scientific Keyboard Input Language (SKIL) to 
Scientific Information Processing Language (SIPL) , and suitable 
character set for CRT typesetting equipment. (RR) 
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Status of Uspekhi Experiment 
July 1969 



Summary 



The main points of the June 1965 Buckland-Weaver proposal to AIP for an ex- 
periment in creating a labeled machine file together with some of the assorted by- 
product publishing outputs, have now been successfully carried out. 

Eighteen pages of journal photocomposition appended hereto were generated 
from a labeled machine file of Uspekhi data by computer programs, as planned. A 
sample of an index page, similarly generated, but from a different journal,* is also in- 



cluded 



These programs are tested and ready to be used in the production of further 

journals and indexes. , 

Costs and procedure breakdowns are shown in Table II-VI. Production details 

will be seen to hav.^ varied in some cases from the 1965 concepts, but not the con- 
cepts themselves. . 

Comparisons between 1965 predicted costs and actual costs m 1969 are shown 



in Tables VII-IX. 

Briefly: 

• The machine file right now costs about 18% more than predicted; it will cost slightly 
less (6%) within the year. 



. The journal page output costs at this point 24% more than predicted, mainly due to 
poorly adapted character generating equipment. CRT character generation promises 
a lowering of these costs significantly-to 27% below. 1 965 predictions. 

• Index output presently costs 14% more than predicted, with the promise of sub- 
stantial savings dependent upon the use of CRT equipment. 

. The systems developed can be applied without changes to immediate Uspekhi pro- 
duction at a cost of $37.9 1 per (English) page. With the installation of on-line PDP-8 
editing equipment (planned for Fall), the page cost drops to $34.00. Further im- 



*The Physical Review. 



provements (provision of an adequate Photon disc and a translation bridge from SKIL 
to SIPL*) would reduce the page cost to $3 1.00. With an investment in a translation 
program and a suitable character set for CRT typesetting equipment, the page price 
falls to $24.00. 

• Assuming the addition to TPS of display-equation programs, SKIL access, and the 
CRT character generator bridge, one can try to predict per-page costs for a primary 
journal on the Uspekhi model: 



Operation 


Per page costs 


(6700 cc) 


Copymark 


$ 1.00 


Keyboard text (4,200 cc) 


6.60* 


Keyboard equations (equiv. 2,500 cc) 


6.00* 


Proofread (1) all data 


4.00 


Proofread (2) storage data only 


.40 


Edit and correct 


2.80 


Edit-display monitor 


1.30 


Generate master file 


1.10 


TPS, including equations 


4.20 


CRT proof run (for author) 


3.30 


Keyboard additional corrections 


2.00 


CRT final run 


3.30 


Page makeup, manual 


3.00 


Proofs, miscellany 


1.00 


TOTAL 


$40.00 


•For simplicity, it is assumed that illustration, tabular material, etc. can be handled 
within the cost limits for a 60-40 text-equation mix. 



Procedures: Creating the Machine File 



Copymarking 

It will be seen from Table II that only $.70 per page has been allocated to copy- 
marking. This is a low figure compared to the line for this item in AIP’s normal pro- 
duction budget. One of the major reasons for this is that the data for Uspekhi has 
already been typeset once in Russia, and a style established during this process. 

It has been found in other experiments, however, that even for primary manu- 
script the amount of copymarking necessary to guide keyboard operators through 
material of this complexity is minimal; certainly far less than in conventional direct- 
setting systems. To a large extent the reason for this is that the need for repetition 

♦See SKIL: Science Keyboard Input Language, V. Weaver, October 1968, a report to AIP. SIPL is 
an acronym for Scientific Information Processing Language, used in Inforonics’ Master File. 
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of complex strings of typographical commands is avoided. The entire command for 
choosing a typeface, leading, measure, and format for a footnote, for example, is fn. 
More than 250 separate information items exist in the Uspekhi program, but the 
identification vocabulary is considerably smaller. Eq, meaning equation, has different 
control meanings depending on whether it occurs in a text or footnote. 

Another reason is that identifier language is chosen for mnemonic value. If Dm 
transfers the keyboarder to the mathematics character set, then a calls out an arrow, 
g calls out greater than, and i calls out infinity. These are standard for all typefaces 
and jobs, and are quickly learned. Since only the general identification is required 
(the specific choice of typeface and point size being made by the computer), there 
would be no point in copymarking such occurrences. 

In practice, about 20 labels (fn, tx, eq, etc.), each with easy mnemonic associ- 
ations, keep recurring, and these can be written in the manuscript margins by people 
withhlittle or no typographical background. 

Codes 

Creation of the machine file requires that all data be labeled as to its nature so 
that future uses, including typographical ones, can be designed to take advantage of 
the labeling. (Author names, for example, can be searched and/or typeset distinc- 
tively, provided they are flagged.) 

In addition, certain purely typographical functions (such as the use of boldface 
for vector expressions, for example) are embedded in the data, since they do not 
follow simple deductive rules. 

It is interesting that in the present case the number of tape codes required to 
label the data on input is almost the same as the number of control codes required 
on the Photon tape (see Table I). They are, of course, quite different. The recogni- 
tion that j3 belongs to the Greek character set and the assignment of a prefix Dg to 
the b code, does not require the same degree of mathematical or typographical 
sophistication needed to recognize that the j3 must be italic because it is embedded 
in a mathematical expression and that a certain Photon matrix position and lens 
must be synchronized to produce the proper point size on output. 

It is surprising that the number of codes needed in the two cases is so nearly 
equal. One objective of input system design should be to make fewer input codes do 
the job. 

The SKIL system (described in reports of March and October 1968 to AIP) 
moves in this direction, cutting the number of keyboard labeling codes by 40%, 
actually somewhat below the number needed to tag the data in machine storage. 

Since the program to translate SKIL input into master-file codes is still missing, 
no consideration of the effect of this potentially powerful tool is given in the cost 
analysis section of this report. 

Keyboarding 

Input languages used in these experiments have been largely independent of the 
kind of equipment employed. Custom built keyboards have been avoided, and input 



systems have been designed to work within the standard 44-keybutton range. Costs 
shown in the tables are based on work with paper-tape perforating keyboards. Ex- 
periments for other purposes have convinced us that magnetic tape keyboards are 
faster and more reliable. Table IV, showing costs for a Fall 1969 production system, 
is built around the use of an IBM Selectric keyboard (#735 I/O), linked to a time- 
shared computer and writing a magnetic tape whose records (typewriter lines) are 
immediately numbered and made accessible to a keyboard-operated editing program. 

This system, known as Astrocomp, is being marketed by Information Control 
Systems, Inc., of Ann Arbor, Michigan. It is scheduled to be installed at VWC in 
August, and we hope it will be operating shortly thereafter. 

Using this system, control language is identified by a key-controlled shift to 
control mode which is accompanied by a ribbon color shift. Otherwise, the key- 
board is the standard IBM Graphics configuration, on 44 buttons. 

Astrocomp is discussed further underlining. 

Listback and Proofreading 

For all input systems in which hardcopy is to be proofread, we have employed 
a listing of the keyboarded tapes to provide an accurate record of codes and to 
identify, via a parity check, as many machine errors as possible. 

The actual method of production for the Martynov article (Table II) included 
a listback of DURA tapes before proofreading. 

The recommended method (Table III) substitutes for this an immediate con- 
version to magnetic tape and a line-printer dump to obtain proofreading copy. This 
step provides a better check of codes, and also a larger character set (and therefore 
fewer characters to be proofread in coded form). 

From a proofreader’s standpoint, an ideal initial printout would involve the full 
character set, and enough formatting to verify that item labels had been correctly 
recorded (including, of course, the assembly of displayed equations when that be- 
comes part of the system). Such a system— involving a preliminary Photon output 
in a typographical form designed exclusively for efficient proofreading, and not in- 
tended for any publishing output— is being used to produce The Physical Review In- 
dex , but has not yet been designed for Uspekhl 

Considering that the computer processing and character-generation time and ex- 
pense for such an ideal step would probably be only very slightly less than that for 
actually producing journal pages, it is doubtful whether this is a practical feature of 
any system designed to handle translation journals. Where a proof of his data must 
be submitted to an author for approval, the costs of such a procedure may seem less 
formidable. 

It will be noted that Table I reflects two complete and independent proof- 
readings of the input data file. If the entire file is to be stored in machine form, this 
is the recommended procedure. If, however, the file be thought of as having two 
parts: one, those items to be stored for later search and processing, and two, those 
items that will be used only once— for typographical production of the primary 
journal; the proofreading technique can be altered. A single proofreading can be 
given the “journal only” portions (which is the present practice under conventional 
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production methods), with a double reading reserved for the “storage” portions of 
the file. This is the procedure shown in Tables III and IV. The assumption is that the 
“journal-only” part of the file will be discarded after output of the journal pages. 

Editing and Master-File Generation 

Tables II and III include the use of Inforonics’ CRT-monitored Edit-Display 
programs for correction of the data. Changes are made as the tapes are fed into the 
Master File. A human editor works from marks previously made on a listback hard- 
copy, and checks his work immediately on a scope display. 

Table IV recommends use of the Astrocomp editing and correction programs. 
These are not unlike the IBM Datatext package (also called the Administrative Ter- 
minal System). Checking is not via CRT, but via computer-controlled printout on 
the Selectric I/O; a slower, but cheaper method. Since the system is not yet in- 
stalled, reliable predictions as to price cannot be made. A small allowance is left 
in the budget for continued monitoring of the read-in operation, with the hope 
that this step will prove unnecessary. 

The Master File format is a standardized coding scheme used within the In- 
foronics Text Processing Service. It forms the basis for all future file manipulation 
and storage functions. 



Cost Tables 



TABLE I. Comparative output character and code counts for the 
Martynov article. 





Total 


Per 

English 

page 


Percent 
of nondata 

codes 


Input keyboard codes, 
Inforonics keyboard 
language 


87,800 


5,100 


18 


Input keyboard codes, 
SKIL 


81,400 


4,750 


12 


Master file (storage) codes 


84,756* 


4,930 


15 


Typesetting tape codes, 
Photon (for journal 
pages) 


87,552* 


5,100 


18 


Output (printing) charac- 
ters and spaces 
(journal pages) 


71,800 


4,200 


— 



•These two totals are computer-counted. Other totals are human estimates. 
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TABLE II. Step-by-step costs for production of the machine file 
for the Martynov article, actual method. 

Per thous- 
and output 





Hours 


Rate 


Total 


characters 


Per page 


Copymark 


1.5 


$ 8. 


$ 12.00 


$ .17 


$ .70 


Keyboard 2 


12.3 


11. 


136.00 


1.90 


7.90 


Listback 


2.4 


4. 


7.60 


.14 


.56 


Proofread (1st) 


6.0 


8. 


48.00 


.67 


2.81 


Proofread (2nd) 


6.0 


8. 


48.00 


.67 


2.81 


Consolidate 












marks 


2.0 


8. 


16.00 


.23 


.94 


Edit display 






102.00 


1.42 


5.93 


Generate 












master file 


.25 


45. 


11.25 


.16 


.66 


TOTALS 








$5.36 


$22.31 



1. A convenient unit for comparison; there is no such thing in the machine file, of 
course. 

2. DURA Mach 10. Paper tape. Inforonics keyboard langu?. 



TABLE III. Step-by-step costs for production of machine Hie, 
recommended and immediately applicable method. 



Per thous- 
and output 





Hours 


Rate 


Total 


characters 


Per page 


Copymark 


1.5 


$ 8. 


$ 12.00 


$ .17 


$ .70 


Keyboard 1 


12.3 


11. 


136.00 


1.90 


7.90 


Convert to mag 
tape, dump on 
line printer 2 


.67 


45. 


30.00 


.42 


1.75 


Proofread (1st) 3 


4.2 


8. 


33.60 


.47 


1.97 


Proofread (2nd) 4 


.5 


8. 


4.00 


.06 


.24 


Edit display 






102.00 


1.42 


5.93 


Generate master 
file 5 


.25 


45. 


11.25 


.16 


.66 


TOTALS 








$4.60 


$19.15 



1. DURA. Inforonics language. 

2. 120-character chain. Partial typographic conversion. 

3. Read full record. 

4. Read storage data only. Storage data consists of information items to be retrieved 
for other purposes than primary journal production. 

5. Full record. 



TABLE IV. Step-by-step costs for production of machine file, 
recommended and applicable in the Fall of 1969. 



Per thous- 
and output 





Hours 


Rate 


Total 


characters 


Per page 


Copymark 


1.5 


$ 8. 


$ 12.00 


$ .17 


$ .70 


Keyboard 1 


11.3 


12. 


136.00 


1.90 


7.90 


Listback 2 


2.6 


5. 


13.00 


.18 


.76 


Proofread (1st) 3 


5.0 


8. 


40.00 


.56 


2.34 


Proofread (2nd) 4 


.6 


8. 


4.80 


.07 


.28 


Edit and correct 5 


2.0 


15. 


30.00 


.42 


1.76 


Edit display 6 






14.50 


.20 


.84 


Generate master 












file 


.25 


45. 


11.25 


.16 


.66 


TOTALS 








$3.66 


$15.24 



1. Inforonics language. Selectric #73 735 I/O. 

2. #', 35 I/O. 

3. Complete record. 

4. Storage data only. 

5. Astrocomp editing programs via keyboard. 

6. Monitored read-in via Edit-Display and scope. 

NOTE: Use of the SKIL system for input was tested earlier on this same data. No 
computer programs are available to read SKIL coding into Inforonics* Master 
File. It can nevertheless be reliably stated that use of SKIL in the production 
system outlined above would result in savings of $1.26 per page for keyboard- 
ing and $.67 per page for proofreading, resulting in a page cost of $13.31. 



TABLE V. Costs for primary journal output from the machine 
file, (Martynov article), actual method. 

Per thou* 





Hours 


Rate 


Total 


characters 


Per page 


Inforonics 
typesetter 
program run 


1.0 


$45. 


$ 45.00 


$ .63 


$ 2.63 1 


Photon 560 run 


5.5 


30. 


165.00 


2.30 


9.65 2 


IBM Composer, 
patches 3 


.65 


11. 


7.15 


.10 


.42 


Check, set or 
steal missing 
characters, 
insert 4 


5.0 


11. 


55.00 


.77 


3.22 


Makeup pages 5 


7.0 


9.50 


66.50 


.93 


3.89 


Proofs, miscella- 
neous 






20.00 


.28 


1.17 


TOTALS 








$5.01 


S20.98 



1. The Infoconics typesetter program now produces paper tape to drive the Photon 
560. Running speed is punch-timited (about 1 10 cps). Conversion to CRT type- 
setting implies magnetic tape, which can be written faster, with consequent redac- 
tion in computer time. Tim would be reflected in a lower price. Since the exact 
amount is not known, the effect is ignored. 

2. No interface program has yet been vritten to driec a CRT c ha racte r -generator from 
the Inforonics typesetter progr am . Such a p rogr am would create access to aenrice- 
buieau composition facilities currently quoting rates per thousand characters less 
than one-fouith thorn for the Photon. It Aonld be borne in mind that none of these 
service bureau operations yet offers a character-set adequate for Uspekhi. The poten- 
tial saving is $7.15 to $7.55 per page. 

3. To translate patches for figure labels and table heads. 

4. The Martynov article was originally chosen to present a difficult ch ar acte r set The 
average Uspekhi article contains approximately one-fifth as many off-disc symbols. 

A price for inserting these characters for the average article would be $ 1 .00 rather 
than $3.22. An investment in a suitable disc would entirely eliminate this item. 

5. This item includes stealing, re-breaking (in some cases), and inserting displayed 
equations from the Russian typeset version. 



TABLE VI. Costs for subject and author index output from the 
machine file, (Martynov article), actual method. 1 



Author and subject sort 


Present 
cost per title 

$ .19 


Mag-tape CRT 
cost per title 

5 .19 


Infoiunics typesetter 


.42 


.42 2 


Photocomposition 


.78 


.30 


Check, dummy, and makeup 


.81 3 


.55 4 


Update master file 5 


.10 


.10 


TOTALS 


$2.30 


$1.56 



1. Taken from experience with The Physical Review Index, now in production. Photon 
560 character generation. 

2. See note 1 to Table V. 

3. This figure subject to reduction with greater work volume; the present load is 16-20 
pages per month, too mall for efficiency. 

4. No consideration a given to the economics inherent in computer-controlled page 
makeup, although (foe indexes) this capabiity is already a part of some systems. 

5. Tins item allows » editor to apply hindsight to subject category loadings prior to 
cumulations^ 



TABLE VD. Predicted vs. actual costs for creation of machine 
file, per page. 





Predicted 

1965 


Actual 

1969 


Percent 

variation 


Actual MarLynov method 
(Table 11) 


$16.25* 


$22.31 


+7 


Recommended method, 
immediate (Table III) 




19.15 


+8 


Recommended, Fall 1969 
(Table IV) 




15.24 


-6 



1. See Jane 1965 proposal, pp. 3 and 4. This figure is the total of items 1-4, 14, 

and 15 ($7.72), multiplied by the Rushm£n^nh page ratio (1.7:1) and ad- 

justed for die Consumer Price Index increase (109.9 to 130.0). Seventy cents 

per page has been added for c o p ymatti ng, which was considered an editorial 

function in 1965. 



TABLE VUI. Predicted vs. actual costs for paged journal-form 
output from machine file, per page. 





Predicted 

1965 


Actual 

1969 


Percent 

variation 


Actual Martynov method 
(Table V) 


S15.10 1 


$20.98 


+39 


Average Uspekhi article 
(Table V, see Note 1) 




18.76 


+2' 


Average Uspekhi article, 
adequate Photon disc 




17.76 


+18 


With CRT character-genera- 
tion (Table V, note 2) 


approx. 


11.00 


-27 



1. See June 1%5 proposal, pp. 3 and 4 and footnote to Table VII. Tim figure com* 
items 5*11 and 13. Note that in 1963 the Linofilm was (a) assumed to have an 
adequate character set, (b) quoted at $20 an hour, versos the $30 an hour now 
being billed for tiie Photon, and (c) predicted to have a throughput rate of 7 cps, 
versus an actual Photon throughput rate di scovered to he 3.5-4 cps. 



TABLE IX. Predicted vs. actual costs for paged index-form 
output from machine file, per article. 





Predicted 


Actual 


Percent 




1965 


1~69 


variation 


Actual Physical Review 


method (Table VI) 


$2.02* 


$2.30 


+14 


With CRT character- 


generation 




1.56 2 


-23 



1. See June 1965 proposal, p. 4 and footnote to Table Vn. 

2. This is a maximum, with no consideration giv en to computer-controtted page make- 
up capacities. See Note 4 to Table VI. 
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INTRODUCTION 

I. In his preface to the Russian translation of 
Robinson and Stokes* monograph Electrolyte Solutions , 
published in 1963, Academician A. N. Frumkin writes, 
“Although the number of experimental and theoretical 
studies on electrolyte solutions at present considerably 
exceeds ten thousand, interest in this field is not waning. 
On the one hand, this is explained by its great practical 
significance both in technology and in laboratory 
practice, and on the other hand, by the difficulties that 
the theory encounters in trying to interpret and generalize 
the experimental facts.** The latter statement is of 
especial interest here, and it can be formulated more 
sharply: the theory of electrolyte solutions is now 
undergoing a definite crisis. As we see it, the reason for 
this is that people try automatically to extend to more 
concentrated solutions the methods that have brilliantly 
proved themselves for dilute systems, while neglecting a 
certain qualitative difference existing between them. 

In fact, as the Debye-Huckel theoryi'i implies, in dilute 
solutions the radius of correlation R { between the 
positions of individual ions (which equals the Debye 
radius = {<0/8 wlcvj*) is considerably greater than 
the mean distance R between particles. Hence, the Debye 
sphere contains a large number of particles at one time, 
and the mean electrostatic-interaction energy per particle 
is small in comparison with 0. In consequence of the 
latter, the spatial distribution of the particles hardly 
differs from that in an ideal gas. Finally, since R c » 
where r, is the diameter of the ions, the contribution of 
the Born repulsive forces to all the macroscopic 
characteristics of the system is negligibly small in 
comparison with that of the electrostatic forces. 

The volume of the Debye sphere - 4 rtg declines 
with increasing concentration as r J2 , while the mean 



•The following notation is adopted here: c is the dielectric constant 
of the solvent, 0 = AT is the temperature e m = k^e is the charge of an 
ion of type a, and k, is its valence, 

v*= J v. 

is the total number of particles per unit volume, and M is the number of 
different types of particles in the system. Hereinafter we shall restrict 
onrselves to treating only binary symmetric electrolytes, for which 

M — 2. k, = —k. = k, and »,=».= ». 



volume per particle v declines as Consequently, e.g., 
in aqueous solutions having 



V ~ 



0.05 — 0.10 



moles /liter 



(I) 



the quantity vgj becomes equal to v. At the same time, a 
substantial rearrangement takes place in the system: 
when v > v^, each ion is now shielded by only one 
counterion. Hence, a pairing of oppositely-charged 
particles seems to occur in solutions of intermediate 
concentrations. In this process, the mean distance /£,_ = 
R. + between ions of opposite type becomes considerably 
less than the mean distance R calculated without account 
of the interaction between the particles.* And although R 
is rather large, as before, in solutions of intermediate 
concentrations. Born repulsion forces begin to play an 
important role because R + . ~ r 0 . Evidently, R+. can 
differ appreciably from R only when the energy of pair 
interaction of the ions >6. 

In going to very concentrated solutions or to molten 
ionic salts, the radius of correlation R c increases again, 
and in the limit it becomes several times r 0 . The structure 
of the melt thus proves to be similar to that of determined 
by Born repulsion forces. 

Thus, the type of distribution of the ions in the system 
undergoes substantial changes with increasing con- 
centration. Correspondingly, the approach in constructing 
a theory of electrolytes of different concentrations must 
also differ. 

II. In 1946, Bogolyubovni was the first to find the 
second (Debye) term of the virial series of a system of 
charged particles by expanding the configuration integral 
of the system in a power series in the small parameter v/vg 
= 3x*. where 



, 8 avkWrl k*e* 

I* — v 

e0 ’ * e0r o * 



x* 4%p, p- 2jivrJ. (II) 



Subsequently Mayer<» and Haga»i have considerably 
refined Bogolyubov*s result by calculating several more 
terms of the series 



p = 1 - IT - y <**>' 2 <2/)H2, -3) -g <* x) * ln <**> 

)=0 



•This effect is manifested more strongly as the maximum value of the 
electrostatic energy. increases in comparison with the energy 0 

of thermal motion. 
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- 11 2 (2, . l)V(2/ — 2) }+°KX*)‘l. CIII) 

)= 2 

where p = /*/ P* = 2vO is the ideal-gas pressure, and 
C = 0.577... is Euler’s constant. One can easily convince 
one’s self by direct substitution that when x > 1* and v 
> vgj (i.e., when x* > 1/3), the virial series (III) 
converges very slowly, or perhaps may diverge completely. 
Hence, condition (I) practically defines the boundary 
between systems for which one can construct a theory by 
expanding the initial quantities in power series in the 
small parameter v/vg (or equivalently, in the concentration 
v), and systems not admitting such an approach. 

At present there are only two groups of methods in 
statistical physics for calculating the configuration 
integral of gaseous and liquid systems: the methods based 
on virial expansions and the so-called methods of 
liquid-state theory, w However, as was shown above, the 
former are unsuitable for constructing a theory of 
electrolyte solutions of intermediate concentrations. 
Hence there remains only one way out: to use for this 
purpose the integral-equations method of liquid-state 
theory .t As will be shown below, one can make some 
progress in this way, mainly because the absolute values 
of the density of the ionic subsystem, which are 
characterized by the dimensionless parameter p = 2xrr 0 \ 
are sufficiently small, as before. On the one hand, the 
latter situation permits us to neglect the nonlinear terms 
in the integral equations, which considerably simplifies 
the entire calculation. Or the other hand, it permits us to 
limit the treatment to a very ample model of the 
electrolyte solution in which the solvent is described by 
its dielectric constant c alone, and the ions are treated as 
hard .^iieres of diameter r 0 Within the framework of 
this model, the energy of pair interaction of particles of 
types a and b is evidently equal to 



<fi™(r) !-0£»(r) t 




0 < r . r 0 , 
r 0 < r oo. 




er 



(IV) 



Both these simplifications lose force in going to 
concentrated solutions or to fused ionic salts. Hence, the 
calculation of the configuration integral for concentrated 
systems is entirely a special problem, which we shall not 
take up in this review. 



•We recall that x a 2-3 for aqueous solutions of univalent 
electrolytes, and for bivalent ones, x a 8-12. 

t Another variant of liquid-state theory, the “free-volume" method, 
which is based on a certain analogy between a liquid and a crystal is 
also unsuitable in this case, since the structure of the ionic subsystem of 
electrolyte solutions of intermediate concentrations is very far from 
crystalline. 




CHAPTER I: THE BINARY DISTRIBUTION FUNC- 
TION OF THE IONIC SUBSYSTEM OF THE 
SOLUTION 

In essence, any integral equation of liquid-state theory 
(and there are a great many of them at present*) can be 
used to describe systems of charged particles. (This is 
because, as a rule, no special restrictions on the nature of 
the decline with distance of the binary-interaction energy 
are introduced in deriving these equations.) Nevertheless, 
all of the progress in the theory of electrolytes has 
involved only several very simple equations^ *-*1 to which 
we shall limit the treatment here. Here we shall begin 
with presenting the theory of Debye and HQckel,«i an 
analysis of which will permit us to reveal the physical 
meaning of the later theories more pictorially. 

1.1 The Dehye-HAckel Theorym 

I. Let us put the coordinate origin at the center of an 
ion of type a bearing the charge e m . Since this ion interacts 
with the other ions in the solution, the mean concentration 
of particles near it will be altered. This will produce a 
spherically-sym metric charge (“Debye atmosphere’’) 
about the central particle, of density 

9a(r)= S e b Xb^ ab (r), (1.1) 

where &(r) is the binary distribution function determining 
the probability of finding an ion of type b at the distance 
r from ion a. If the two particles are so far apart that they 
no longer interact, then all correlations between them 
vanish, and 9+ becomes unity, while q£r) becomes the 
constant 

fl.(°°)a S e b v b' 

Obviously the latter quantity must equal zero, since 
otherwise an electric field of non-zero intensity would 
exist at great distances from our selected particle. This 
would give rise to an electric current. Hence, the 
condition of neutrality 



21 ( 1 . 2 ) 

must be fulfilled in order that the system of charged 
particles be at equilibrium. 

We can derive another, very important formulation ot 
the condition of neutrality from the well-known formula 
of electrostatics!' 0 ) 



<Pa (r) = 





dV 

*|r-r'| ’ 



(1.3) 



which describes in this case the potential distribution <p t (r) 



•The most complete review of the results obtained in the theory of 
ordinary liquids is given in 0*1. 
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near the central spherical particle. By using Green’s 
formula to transform the integral appearing in this 
expression, and taking into account the fact that q.{ «> ) = 
0, we obtain 

r t 

m 

+-^-{*« + 4 ji J g«(r)r*dr} . (1.4) 

Now we note that any macroscopic system can have finite 
values of its thermodynamic parameters only when 






Correspondingly, as is known, q . must approach zero 
faster than r-\ and tp. faster than r->. 

However, the latter is possible only when the expression 
enclosed in curly brackets in (1.4) is zero. That is 

m 

e a \- 4n \ e b v b Vab (r) r* dr 0. ( 1 .5) 

In other words, in order that the binary distribution 
function decline rapidly enough toward infinity, the sum 
of charges of the central ion and its atmosphere must be 
zero. Fulfillment of condition (1.5) provides for shielding 
of the field of the central ion, since only in this case does 
the integral term in (1.3), which represents the potential 
producing the Debyr atmosphere, “throw out” a term 
“cancelling” the potential ejtr of the central ion. As will 
be shown below, such a situation is typical of the integral 
equations describing systems of charged particles. 

II. The Debye-Huckel theory is based on the idea of 
using formula (1.3)t to calculate Sjir). Evidently, to do 
this we must assume that = .V «*(ip.(r)), since otherwise 
(1.3) does not form a closed equation. The concrete form 
of the sought relation can be found as follows. 

First we note that the definition of the binary 
distribution function implies that it can always be 
represented in the formal V^(r) = y(r)G+(r), where the 
coefficient y (r) = exp [— ipw(r)/©] does not depend on 
the charges of the ions. Hence, .V* and <P. can be related 
only by way of the function G*. Second, S + must by 
definition be a symmetric function of the charges e. and e k 
of the particles, w On the other hand, we can easily derive 
from (1.3) and the neutrality condition (1.5) the fact that 
~ e,. This implies that G* can depend only on the 
product e k ip m . Third, since G* is dimensionless, it can be a 
function only of the dimensionless potential etf./® = 
(ejeJ t®r 0 )\p{r) = (e«e*/fcV0xW- 

Thus, if is determined by the value of <p„ then this 
relation must have the form '§+ = yGJ^ejeJ Mx^l- Let 



•For example, the density fluctuations, which are proportional to 
Jj* — IjrVr, go to infinity when this condition is violated. 

tOr, equivalently, Poisson’s equation, of which (1.3) is the solution. 




us expand the unknown function G^ in a series, and limit 
it to the first two terms. In this approximation. 



Sab (r) ' - Y (0 [ * - jQ- *( r ) j • < 1 1 - 6 > 

The unknown expansion coefficient B 1 appearing here is 
in general a function of the dimensionless parameters x 
and k of the problem. Debye and Huckel assumed it to be 
unity. As will be shown below, the more exact theories 
confirm this assumption. 

III. Taking £* = 1, let us substitute (1.6) into (1.3), 
taking (1.1) and (1.2) into account. Consequently, (1.3) 
acquires the form 




where t = r/r 0 and y(/) = 0 when 0 ^ / < 1;* or 1 when 
1 S / ^ oo . Let us calculate the integral appearing in 
(1.7) using Green’s formula, using the condition of 
neutrality (1.5). Thereupon, (1.7) is reduced to 



OO JO 

hx* «• <■ *> 



Now assuming that ^ we transform the 

integral equation (1.8) into a differential equation: 



tn." (/) — x*m ( l ) = (l. (1 .9) 



It directly follows that m = Ae~ u . Finally, we determine 
the constant A from the condition of neutrality ( 1 .5), and 
we obtain 

&W“v(0{i-^r*nr}’ (U0) 

Upon substitution into the expressionf'i 

U (F, 0)-f/ o (<->) 

r I? V ?^<f'*b{r)S ab (r)'i3ir'-dr. (1.11) 

0 1 <a, 6cJM 

which relates the internal energy U of the system to S+, it 
gives ($* being defined by formula (IV)) 

U{\\ ») ■ (1.12) 

This is the final result of the Debye-Huckel theory. 

IV. We shall make some remarks on the formulas 
derived above. 

First, the Debye-HQckel theory is based on the 
hypothesis that, in the region r > r 0 , the distribution 
function arises only from the Coulombic interaction 
among the particles (since it is assumed that G* = 
GJetf./G). Evidently, the latter implies complete neglect 



•Of course, the latter is true only of systems of hard spheres. 
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of short-range forces at distances exceeding the diameter 
r 0 of the ions. However, as is known, ^ 1 for r > r 0 , 
even in a system of uncharged hard spheres. Here the 
thickness of the spherical shell r 0 g r g R r W in which the 
positions of the individual particles are correlated 
increases rapidly with increasing concentration. »•"> 
Conversely, the thickness of the spherical shell Rw — r 0 
manifesting electrostatic correlation declines with 
increasing concentration. The decline approximately 
follows the law Rw — r 0 ^ r 0 /n ~ 1 Vv, as (1.10) 
implies. Hence, the condition /? r «« — r 0 <<c /?,<*» — r 0 , 
which has to be fulfilled in order that we can take G to 
depend only on the electrostatic potential holds only 
over the limited concentration range p < p mn . If we use 
the more precise values of the radii of correlation R c li) and 
Rw given in i 912 ’, we find that for aqueous solutions of 
univalent electrolytes for which 2.0 ^ x = 3.0 (see 
below), the value of p m « 0.2. This corresponds to 
concentrations of the order of 2-3 moles/liter.* 

Second, the Debye-Huckel theory is based on 
expanding the function G+ = G^ejeJk^xty J in a 
power series in x*lAr), and dropping terms of the order of 
(x^) 2 and higher. Evidently, the latter is possible only 
when \4Ar) « 1 throughout the region r ^ r 0 . Hence, 
the results obtained become more accurate for smaller 
other conditions remaining the same. And since x * s 
proportional to k 2 , it follows that the Debye-Huckel 
theory is inapplicable to solutions of intermediate 
concentrations for bivalent electrolytes, for which x = 
8-12, as is known.* At the same time, the Debye-Huckel 
theory should give quite satisfactory results for univalent 
electrolytes, for which the product x^ is less than unity 
for almost all r > r 0 .+ 

Third, the Debye-Huckel theory has the fundamental 
defect that it doesn’t permit one to obtain full 
information on the distribution function. Indeed, we note 
that, since the charges of the particles enter into the 
configurational energy U N of the system only in the form 
of pair products, therefore G* = GJ,e&). Hence, let us 
represent G+ in the form 

G a l<~' I — " ga b^'-'a^h) 1 \gab i^a^b) fnt ( ^a^b) I 

■f Y (Vb) gab ( — fa*^)], (1.13) 

where we have written the term unity separately in order 
to emphasize the fact that G+ — 1 when r — * ». 
Evidently, for a symmetrical binary electrolyte, for which 



•It is usually assumed that the Debye-Hflckel theory is applicable 
only for v ^ 0.01 moles/liter. This statement is based on a 
misunderstanding, since it results from comparing formula (1.12) with 
virial series such as (III). However, such a comparison cannot serve as a 
satisfactory criterion, since when v > v^, the virial series generally 
diverge, while (1 .12) continues to hold true. Insofar as we know, no one 
has yet made any special estimates of the limits of applicability of 
formulas (1. 10) and (1.12). 

tOne can easily convince one’s self of this by considering Eq. (1.1). 

[14 



o 




i, = —e. = ke, the expression within the left-hand 
square brackets changes sign as we go from G ++ = (7_ _ to 
G + _ = G + , while the second square bracket remains 
unchanged. This implies that the binary distribution 
function can always be written as 

:U V (0 (l --£ixMV) : «(')} . (1-14) 

where we have written the coefficient ( — x) separately 
only for convenience. Evidently, since 

to fKS^-D-o. 

f-too 



liin t*M (/) - - 1 im / 3 V? (() 0. (1.15) 

t ~+ f X> t — ♦oo 



Substituting (1.14) into (1.3) and (1.5), we obtain 

*<'>-■ r + lfS V(OAf(<') r j?v T = 0, (1.16) 

TO 

1 - X 2 \ M 0. (1.17) 



Thus we see that the function R(t) completely vanishes 
from all the electrostatic equations, owing to the 
condition of neutrality (1.2). Nevert* .»ess, if we could 
find from (1.1 5)— ( 1.17) the exact value of the function 
M(t), then this would suffice for correct calculation of the 
internal energy U of the system in the case in which the 
pair potential $*(r) is given by relation (IV). In fact, by 
substituting (IV) and (1.14) into (1.1 1) we obtain 

oo 

U(V. 0) Uo(Q)— M(t)tdu (1.18) 

i 

That is, U(V ,©) does not depend on R(t) at all. At the 
same time, the pressure of the system 



/'O'. «)- J*id(r, 



»)— i ^ 

0 1 'a. fcC-W 



(Wqb ( r ) 

dr 



- 'bah (r) 4nr 2 dr (1.19) 

cannot be found by using only one function M , since 
after substituting (IV) and (1.14) into (1.19), we get the 
expression 

oo 

P -£=1 \ + «(l)l, (1-20) 

from which R(t) does not drop out. Hence, (1.10) gives 
V(V, 0) correctly, but P{V, ®) incorrectly. 



1.2. The Kirkwood-Poirier Theory w 

I. Following Debye and HQckel, we set \p(t) = M(t) in 
(1.16). Here the physical meaning of Eq. (1.16) is changed 
considerably. Indeed, it formerly gave the electrostatic 
potential rpit) at some arbitrary point t not directly 
associated with any particle of the system. However, after 
\f/ has been identified with M , it now gives the probability 
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of finding the particle b at the distance t from the 
“central” particle a (since the correlation function M is 
part of the binary distribution function y ab ). Hence, while 
the integral found in (1.16) formerly described the field 
created by the N — 1 particles at the point t (i.e., all the 
particles of the system but the central one), now it must 
describe the influence of N — 2 particles on the pair of 
particles a and b that we have chosen, whose positions 
are fixed. Correspondingly, we must perform the 
integration in ( 1 . 16) over the entire space V, excepting the 
volume occupied by these particles. However, the factor 
y(t') within the integral already takes into account the 
excluded volume occupied by the particle a. Hence, we 
have yet only to take into account the volume of the 
particle b. We can do this by introducing the factor y ft - 
t') into the integral. Consequently, (1.16) is converted 
into the equation 

M (/) - j + £ \ M «') V «') V (I l - *' I) rpTj ’ <*. 

* (1.21) 

which is the basis of the Kirkwood-Poirier theory. We see 
from the arguments given above that it is a natural 
refinement of the initial equation (1.7) of the Debye-Huckel 
theory. 

II. The integral over the volume occurring in ( 1 .21 ) can 
be transformed by Green’s formula into integrals over 
the surface of spheres having centers at t = 0 and at t. 
Upon substituting M = m" t, the latter are easily 
calculated. Consequently, when t > 2, Eq. (1.21) 
becomes 

00 

2<f<00, — — y—— J* ^ 1 — X 2 \ A/ (/)/*<//} 

1 

— (m(/ 1) m(t 1)] 0. ( 1 .22) 

If we let r — ► ® here, and take (1.15) into account, we 
come again to the condition of neutrality (1.17), which 
must be satisfied to remove the l/y xxt term from Eq. 
(1.21). After simplifying (1.22) in /, we finally obtain (cf. 
(19)): 

2 t .oo, m n (t)-^-[m(( -■ 1) ;>/*(/ — 1)1 0. (1.23) 

If the distance between the particles t < 2, then the 
spheres interact (Fig. 1). Correspondingly, the integral 
must be taken only over the outer surface of the figure 
formed. The end result of this is the expression 

\<t< 2, m" (/) - y x 2 m (/ - 1 ) j v.-m ( 1 ) 

- x* ( 1 — m' (1). (1.24) 

The region 1 ^ ^ 2 is marked out physically because 

there is always a “forbidden” zone between two particles 
of diameter r 0 when r 0 ^ r ^ 2r 0 . The third particle over 
whose coordinate the integration in (1.21) is performed 
cannot enter this zone (see Fig. 1 ). 




FIG. 1. Diagram of the collision of three particles with r 0 g r ,2 g 
2ro- 

We can easily convince ourselves by direct substitution 
that a special solution of Eq. (1.23) is e-«, where a is the 
root of the transcendental equation 

a 2 - x 2 ch a. (125) 

Since (1.25) has an infinite number of roots, the general 
solution of (1.23) can be written as 

oc 

/«(/) V^-V (1.26) 

» I It" 

The unknown constants A , occurring here must be 
determined from Eq. (1.24). Upon substituting (1.26) into 
the latter, it gives 

oo 

1 / 2. ^ .-I -- ~ x- {m (1 ) (/ — 2)m'(1)J. 

. i “ (1.27) 

where A* = (AJA^fa? - V 2 K 2 e«i]. In essence, Eq. (1.27) 
is an expansion of the known function (within the 
accuracy of the constants w(l) and m\ 1)) on the 
right-hand side in a series of the characteristic functions 
e-a^ of Eq. (1.23). In principle we can find from (1.27) all 
the A , but A „ which actually drops out of ( 1 .27), owing to 
homogeneity. The constant A, can be determined from 
the condition of neutrality (1.17), which gives the 
following when (1.26) has been substituted into it: 

A , (x 2 (m (i) — m' (I))! -1 . 0-28) 

The study of the transcendental equation (1.25) made 
in "«i showed that it has only two real roots a u = a , y, the 
remaining roots are complex conjugates: a k = a k ± iu j*. 
When k < 0.5, a, % «, a 2 ^ 5.5, and the values of the 
other a k are even greater. Hence, when k £ 0.5, we can 
drop all the terms of the series in (1.26) but the first, to a 
sufficient degree of accuracy. As a result we obtain the 
Debye formulas (1.10) and (1.12) for .V* and U. 

III. We shall now make some remarks on the 
Kirkwood-Poirier theory. 

First, we stress the fact that its starting equation (1.21) 
was initially derived directly from the Gibbs canonical 
distribution, rather than from the electrostatic equation 
(1.3) with the substitution of (1.6), as was done above. 
Hence, the problem of the value of the expansion 
constant B 1 in (1.6) does not arise at all in this case. We 
see from the . derivation given in w that B 1 = 1 . 
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Second, in Eq. (1.10), which is also a consequence of 
the Kirkwood- Poirier theory, let us transform to the limit 
of uncharged hard spheres having e a = e b = 0. Then 
proves to equal y(r). Upon substitution into (1.19), this 
gives the expression p = 1 + %p, which amounts to the 
first two terms of the virial series. i"i We can easily 
convince ourselves by comparing p = 1 -I- with the 
more precise equations of state that this expression gives 
satisfactory accuracy only when p < 0.2. For r 0 ^ 3-4 A, 
this corresponds to a concentration v s* 2-3 moles/liter. 
As we should expect, this value coincides with the limit of 
applicability of the Debye-Huckel theory (see above). 

Third, Kirkwood and Poirier derived Eq. (1.21) by 
expanding the initial quantities in a power series in x. and 
then dropping the terms proportional to x*» where n ^ 
2.* This is the very reason why (1.10) gives correctly only 
the first two terms of the virial series of the system of 
charged spheres (III) depending linearly on x* the 
remaining terms of the series (III) depending on the 
higher powers of x cannot be found by using (1.21). 

Fourth, Eq. (1.21), just like (1.7), determines only the 
“electrostatic” component M(t) of the binary distribution 
function; its “short-range” component R(t) cannot be 
found within the framework of the Kirkwood-Poirier 
theory. 

13. A Theory Based on Bogolyubov’s Equations!**) 

I. Let us apply the operator V to (1.16), and assume 
that ^ = M, following Debye and Huckel. Then by 
following the same argument as in deriving the 
Kirkwood-Poirier equation (1.21), we obtain 

T.l/ (t)-V (1) 

f ■£ \ w (O y (O y <1 * - »' ;> ▼ f rt 4rr ) ~ °- 

V (1.29) 

This expression is the equation of balance of electrostatic 
forces acting on the particle b .t We can easily convince 
ourselves of this by noting that in ( 1 .29) the term V ' is 
proportional to the electric field intensity produced by 
the central ion a at the point t, while the integral is the 
mean electric field intensity produced by the remaining N 
— 2 particles of the system at the same point. Evidently, 
the force balance written in this way cannot under any 
conditions be considered complete. This is because it 
does not take into account the non-Coulombic forces 
causing the impenetrability of the particles. Hence, a 
natural generalization of the theories developed above is 
to go over to Bogolyubov’s equations 

4- $ ab V x 0 ab + \ 2 v c V abc V,</>«cd 3 r 3 ---- 0, 

v i<ccm (1 .30) 

*Ste also the original derivation of (1.21) given by Falkenhagen and 
Kelb^.nsi 

t Or, equivalently, the forces acting on the particle a, since the latter 
arc equal except for sign to those acting on the particle b. 



which is the equation of balance of all forces. In fact, let 
us divide each term in ( 1 .30) by In the expression that 
we get thus, the term — V,4>^ will represent the total force 
exerted on the particle a by the particle b , whose position 
is fixed, while 

-IV v c($abclSab) V i^acd 3 r 3 



is the total average force exerted on a by the remaining N 
— 2 particles of the system (here $ aht is the ternary 
distribution function). Correspondingly, the term — V,{® 
In & will represent the “force” of the thermal 
movement of the particle a balancing the potential forces 
acting on it. We see from the derivation of Eq. (1.30) 
given in w that it follows strictly from Gibbs* canonical 
distribution. This is the very reason why we can write it in 
the form of a condition of constancy of chemical 
potential of the particle a in the system in which the 
position of particle b is fixed. 

II. Among the infinite multitude of solutions of 
Bogolyubov’s equation (1.30), the only ones that fit the 
canonical distribution (and hence have physical meaning) 
are those that satisfy the conditions of normalization, 
symmetry, etc.' 2 ’ >» In essence, the latter play the role here 
of ordinary boundary conditions as used in solving 
differential equations (since both have as their purpose 
the selection of a given solution from the entire class of 
possible solutions). Hence, it is natural from the outset to 
try to put (1 .30) into such a form that the imposition of 
additional conditions (or at least part of them) would be 
superfluous. 

The definition of the distribution function implies! 1 u 
that and IS abt can always be represented in the form 
&Jr t2 ) = y(r l2 )[l +gj, .W(r, 2 , r 13 , r 23 ) = y(r l2 )y(r i} )y(r 2} ) X 
[1 +Gw].* Here the unity terms are written separately to 
emphasize the fact that 



and 



lim 






lim 



r i2 ,r i3 ,r 23 ' 



if abc I ■ 



The condition that the correlation should decline also 
implies that if, in the group of three particles <z„ b 2 , and c 3 , 
one of them (e.g., c } ) moves infinitely far away, then IS ahc 
becomes equal to S ab . Hence, when r, 3 , r 23 — * ® , the 
function G ahc — » g ah . Taking into account the symmetry of 
with respect to permutation of particles, we obtain 



•*Ai6 - Y 1 2 1 1 : 



Vab< Y 12Y 13Y23 1 1 -gilb gne ‘ gb': tfabcli (1.31) 

Evidently, the correlation function g tbc describes here the 



•Here as before, y — exp [— cp <,, /01- 



[16 



STATISTICAL THEORY OF ELECTROLYTE SOLUTIONS 



177 



non-linear effects that arise in the simultaneous 
interaction of all three particles. Substituting (1.31) into 
(1.30), we get a system of two equations for the functions 

g ++ = g__ and g + _ = g_ + : 

*Wigab+(l+gab) V l^ ,) 

1 % 1 

+ \ 2. 1 g f gbc -r S*6cl 

XYi3Y2 3 V|0 ac dV 3 -O. (1.32) 

This way of writing Bogolyubov’s equation (1.30) has the 
advantages that it ensures satisfaction of the conditions 
of normalization and symmetry, the condition of 
declining correlation, and the condition that the 
distribution function be bounded at the zero-point,* for 
any approximation of the unknown function g^. 

III. Now we recall that each ion is shielded by only one 
counterion in electrolyte solutions of intermediate 
concentrations. Here the mean distances between the 
thus-formed particle pairs are rather large, as before. 
This means that in these systems the simultaneous 
collision of three or more particles is a rather rare event. 
And if this is so, then we can completely drop the term 
git* in (1.32) describing the non-linear effects in ternary 
collisions. Using (1.14) to transform from g. b to the 
functions Af(/) and /?(/), we obtain instead of (1.32) the 
system of two equations: 

V.V(/)-v(|) 

£ \ .1/ (Oy<Oy(|»- »' I ) v r ,— rf dV 

Y<0*v(|i“i'i>rfv /'(')* (4) . 

(1.33) 

'■«<'•>— C- i n • no : «(') 

t 

V' 

rt (I « — t' |) Y (' ') (!* — *' I> «(’/' - X ! { -»(')* (4) 

+ ~\ fi)Y(|t-f|) Y (OT--! r , 7 rfV. (1.34) 

The left-hand sides of (1.33) and (1.34) are the equations 
of balance of the electrostatic (cf. (1.29)) and short-range 
forces, respectively. The right-hand sides describe the 
interaction effects. 

IV. We shall restrict ourselves in solving the system 
(1.33) and (1.34) to the case of not very great 
concentrations, for which p « 1. Dropping all the 
terms proportional to p in (1.33) and (1.34), we obtain 



•The latter condition is ensured by writing the factors y separately in 
and 



*M (/, — V (4) 

T 14 ? M (/') Y (O Y (I I ~ f I) ▼ nVi <Pf 

V 

(1.35) 

oo 

T«(0-, x W(0ri, 1/ (!) - x* \ -4P-./T.0.36) 

By assuming in (1.35) that M = m^/)//, and 
transforming the integral contained in it by Green's 
formula, we can reduce the system (1.35) and (1.36) to a 
single equation: 

— 4r( m '( , + 1 >-m'(<-l)l-0. (1.37) 

In going from (1.35) and (1.36) to (1.37), we have taken 
into account the condition of neutrality ( 1 . 1 7). As before, 
the latter must be fulfilled to remove from (1.35) the term 
V 1 ^, which declines too slowly toward infinity . When t 1 > 
*X. we can drop the term containing x 2 in (1.37).* As a 
result, we obtain the following equation (cf. (1.9) and 
(1.23)): 

t* = y xx < t~< oo, 

m IV (0— ~ [m ' (t -f 1) -m' (/ - 1)1 - 0, (1.38) 

Evidently, the general solution of the latter has the form 

oo 

m (t) - V As-*' A t m{l), a\ - x. 2 ^ . (1.39) 

»=i n ' 1 

The constants A , , z ^ 2 appearing here must be 
determined from the condition that Eq. (1.37) should 
vanish identically in the region 2 £/&/*, and the 
equation obtained from (1.35) and (1.36) should vanish in 
the region 1 ^ ^ 2, which is considered separately, as 

before. Just as in the Kirkwood- Poriier iricr theory, the 
constant A, is found from the condition of neutrality 
(1.17), which leads to the expression (1.28). 

The study of the transcendental equation (1.39) made 
in i*i shows that when x < l,t the first root of this 
equation a, ^ k, while the real part of the remaining 
roots a, £ 6, / ^ 2. Hence, when k £ 1 , we can drop all 
the terms of the series (1.39) but the first. Consequently, 
the expression for the binary distribution function takes 
on the form 



•This estimate can be obtained by substituting the special solution 
e- fl i*of Eq. (1.38) into (1.37), and matching the terms in * 2 and k 1 . 
t Rather than k < 0.5, as in ihe Kirkwood-Poirier theory. 
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ou 

"■ ' \— ,-*} • (140) 

/ 

In distinction from the Debye formula (1.10), it is 
positive for all values of t and x- However, upon 
substituting ( 1 .40) into the formula ( 1 . 1 1 ) for the internal 
energy, we get the Debye expression (1.12) as before. 

V. We shall now make some remarks on the theory 
developed here. 

First, if we go to the limit e. = e b = 0 in the initial 
system (1.34), (1.35), then we obtain 

V/ -'(0--£ \ 11 : //('H H(t') 

V 

r ^ (I t - 1' 1)1 Y (/') V Y (I t - 1' I) d 3 /' ---0. (1.41) 

It was shown in «'« that (1.41) gives correctly the first 
three virial coefficients, but gives the fourth one with a 
small error, which amounts to only 3.5 percent for a 
system of hard spheres.* The more precise solution 
obtained in i |2 i showed that (1.41) gives quite satisfactory 
results up to p « 0.6- 1.0. For r 0 = 4 A, this corresponds 
to a concentration of 20 moles/liter or more 1. The latter 
overlaps the entire region of electrolyte solutions of 
intermediate concentrations with room to spare. 
However, we cannot say the same for the system (1.35) 
and (1.36), since it was derived from (1.33) and (1.34) by 
dropping all the terms of the order of p. In this 
approximation, (1.41) correctly describes only the first 
two virial coefficients. Correspondingly, the region of 
applicability of the solutions obtained cannot exceed 2-3 
moles/liter. 

Second, as the derivation given above implies, the 
fundamental system of equations (1.33) and (1.34) 
involves no expansions whatever in the parameter x- This 
is the very reason why it gives correctly the first three 
virial coefficients of the system of charged hard spheres. 
We can convince ourselves of this as follows. In the case 
of dilute enough solutions at small distances, we can 
neglect in (1.32) not only the non-linear effects in the 
collective interactions (i.e., the term g mlK ), but also 
generally all the collective effects described in (1.33) by 
the integral term (since the latter are small in comparison 
with the pair interaction described by the term V^. t <«»). 
Thus we directly find that at small distances 1 + = 

exp [( — e,ej t®)(\ /r)\. Correspondingly, for great 
distances with k * « 1, we find from (1.40) that 1 + 

= 1 — (e.ej t © j • { [exp (— r/r^)]/r j. Combining the two 
results together into a single interpolation formula, we 
arrive at the well-known Tyablikov-Tolmachev» 7 >*i 



expression 



$ab (r) -- y ( r ) exp [ 



ee r J • 



(1.42) 



Upon substituting (1.42) into (1.19) and calculating the 
integral by expanding the exponential in a power series in 
(^*/f©)-[exp (—r/rg)/r], we obtain p = 1 — (x*/6) — 
‘/2(x*) 2 X 2[x^-y(2/)(2y — 3)] + 0 (k 3 ). This exactly 
corresponds with the first three terms of the series (III). 
Since here the k 2 term contains all the powers of x* 
evidently dropping the g alK term in (1.32) imposes no 
limitations on the value of the parameter x- 

Since the formula ( 1 .42) can also be derived from the 
“abbreviated” system (1.35), (1.36), the latter remark 
fully pertains to it as well. However, this does not mean 
that the expression (1.40) for the binary distribution 
function found by solving the “abbreviated” system is 
valid for all values of x* As was shown in pi, the result of 
dropping all the terms of the series of (1.39) but the first 
is that (1.40) describes only aqueous solutions of 
univalent electrolytes, for which x does not exceed 2-3; it 
is not valid for multivalent electrolytes. 

Third, as we see from (1.40), the equation of balance 
of all forces (1.32) completely determines the binary 
distribution function. This permits us directly to calculate 
not only the internal energy U of the system, but also the 
equation of state P(V, ©). When ^ and arc given by 
the formulas (IV) and (1.40), the latter has the form 



P 1 



i x* 

(i 1 -| x 



2 






If we know P and U, we can test the theory for 
self-consistency. In fact, we know from thermodynamics 
that P and U are connected by the relation n»> 



(a(//aF)© = 02(a/a0XP/ev. 



Hence, 

P = ©| J%U/d V)(dQ/Q 2 ) + const | .* 

If we substitute into this V from (1.11), we obtain 

P- l- *>-(1 ; x) 1 — 21.1(1 i x)| r|f>. 

(1.44) 

This expression differs greatly from (1.43) in external 
form. However, as we see from the data of Table I, which 
gives the values of A = — pc •*» calculated for x = 2, 

the quantitative difference between the two formulas 
practically nowhere exceeds the experimental error d = 
±1 X 10- 3 . The fact that A nevertheless is not strictly 



*We recall that (1.7) and (1.21) give correctly only the first two virial 
coefficients. 

[18 



•The constant appearing here m ist be assumed equal to Vip for the 
case of hard spheres in order to take into account the excluded volume 
of the particles within an accuracy of terms of the order of p, inclusive. 
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zero is most probably to be explained by the fact that 
(1.40) is an approximate solution of the system 
(1.35), (1.36). Thus, the initial equation unambiguously 
determines the thermodynamic parameters of the system 
of charged si/heres* This is quite natural, since the 
conditions of normalization, correlation, etc., in turn 
determine unambiguously the linear terms of Eq. (1.32). 

We see by comparing the equations for the distribution 
functions (1.7), (1.21) with (1.33), (1.34) that each of the 
latter is a natural refinement of the former. Without 
vitiating the fundamental corollary of Eq. (1.7) (the 
Debye-Hiickel formula (1.12) for the internal energy), 
these refinements nevertheless show that the latter is valid 
only for k < 1 and x ~ 3. A direct calculation of the value 
of U(V , ©) by the Monte Carlo method confirms this 
conclusion. < 3S i 

CHAPTER II: THERMODYNAMIC FUNCTIONS OF 

SOLUTIONS AND COMPARISON WITH EX- 
PERIMENT 

The formula (1.12) for the internal energy of the ionic 
subsystem of the solution was first derived by Debye and 
Huckel as early as 1923.m In the 40 years that have 
passed, all its consequences have been subjected to 
thorough experimental testing. It has turned out that the 
theory does not satisfactorily agree with experiment at 
any appreciable concentrations unless one introduces 
empirical corrections . bmu Up to now, the observed 
disagreement has usually been explained by saying that 
the theory itself is crudely approximate and semiempirical. 
Although this viewpoint has become widespread now, it 
is wrong, since as was shown above, Eq. (1.12) is a direct 
consequence of the Gibbs canonical distribution. But 
then the question arises anew: what then is responsible 
for the stated discrepancy? 

To answer it, we recall that the theory involves the 
dielectric constant c of the solvent, which arises from 
averaging the configuration integral of the system over all 
possible states of the solvent molecules. Evidently, 
the result of this averaging cannot fail to depend on the 
temperature © of the system and the concentration v of 



•One can show that Eq. (1.41), which is a special case of 
(1.33), (1.34), also unambiguously determines the thermodynamic 
parameters in the case of a system of uncharged particles. 



the solute (and of course, on the other parameters of the 
electrolyte and the solvent, such as the valence of the 
ions, the diameter of the particles, etc.). Hence we can 
state a priori that always c = «(©, v). 

It was noted long ago that one must take into account 
the temperature-dependence e = «(©) to explain certain 
properties of extremely dilute aqueous solutions of 
electrolytes.* On the other hand, the concentration- 
dependence i(v) was said to play no role. However, such 
a disparity between @ and v only stems from the fact that 
the concentration parameter c, = (p/c)(9c/9p) is always 
small in dilute solutions, while the temperature parameter 
«© = (®/«)(9«/9©) for water proved fortuitously to be of 
the order of unity (more exactly, cq = — 1.4.i»i). The 
situation alters with increasing concentration, since then 
the value of c, also rises. For aqueous solutions it attains 
a value of ^ 0.2 at v = 1 mole/liter. i J5 i If in addition we 
take into consideration the fact that c, occurs in the final 
expressions in the form of the product x«»t and that 
always x £ 2-3, then it becomes obvious that we must in 
no way neglect the concentration-dependence c = c(p), 
even at v > 0.1 -0.2 moles/liter. Unfortunately, this 
feature hasn't received due attention in calculating the 
thermodynamic characteristics of electrolyte solutions of 
intermediate concentrations. On the one hand, this has 
given rise to some inner contradictions in the ther- 
modynamics itself (see below). On the other hand, it has 
made it impossible to attain good agreement of theory 
with experiment. 

2.1. Thermodynamics of Electrolyte Solutions 

I. Averaging the configuration integral of an electrolyte 
solution over all possible positions of the solvent 
molecules is equivalent to going from the real solution to 
a model system consisting of N charged particles, each of 
which moves freely (i.e., frictionlessly) in a continuous 
medium of dielectric constant c. However, the previous 
chapter essentially treated not this model, but a gas of 
particles bearing the effective charge el V7, and moving in 
a vacuum. Hence we must first of all analyze to what 
extent the results obtained above can be used to construct 
a thermodynamics of electrolyte solutions. 

All calculations of the binary distribution function of 
an ion gas are based on the assumption that its 
configurational energy U N equals the sum of the pair 
potentials 4>,*(r) given by relation (IV). In order that U s 
should have the same form for the model system as well, 
we must assume: a) that the charge within the ion-spheres 
is distributed with spherical symmetry, and b) that the 
dielectric constant c<o of the medium within the ions 



•Without this, for example, one fails to explain not only the 
concentration-dependence of the heats of dilution, but even the sign of 
the effect itself. > 24 > 

tSee, e.g., the expression (2.7) for the free energy, which contains 
the term !6(x - *o) a V$xo«,. 
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equals the dielectric constant < of the solvent.* If both 
these conditions are satisfied, then we can assume that 
the above-derived formulas for 27 * are equally valid both 
for the ion gas and for the model. However, the 
expressions for the thermodynamic functions of the ion 
gas derived by using (1.1 1) and (1.19) cannot be extended 
to electrolyte solutions. This is because, generally 
speav.ng, these formulas do not define the relation 
between the binary distribution and the macroscopic 
characteristics of the model that we are studying. The 
latter can be seen especially well from the example of the 
formula (1.19) for the pressure, which implies that the 
solvent influences the equation of state of the solution 
only by way of the ionic subsystem. However, as is 
physically evident, the contribution of the ionic 
subsystem to P is generally very small, and fundamentally 
P depends precisely on the properties of the solvent. 

In order to establish the relation between the binary 
distribution function and the macroscopic parameters of 
electrolyte solutions, we shall start with the well-known 
expression for the internal energy of the system 



\ Of 



r , e ot\ ■ ,h‘ _ 

V 1 - ) T: 

V 



Since in our case 



( 2 . 1 ) 



E= Ee,, 

where E, = E,(lr — r,|) is the field of the rth ion, then 



e 2 = Ee, 2 +Ee / e / . 

i i=i 



Substituting this expression into (2.1), we obtain the 
following for the electrostatic component of the 
configurational energy of the solution: 



lT 



A ( i . / V V — ) 

- \ ' e l — . ff 0 ' — t'ij ) ’ 



I • : i , v 



i<i, j<y 



( 2 . 2 ) 



•If «« * «. then the solution of Poisson's equation = 

. . . r v ) for the system of charted spheres, where the r, are 
the coordinates of the rth ion. cannot be represented in the form 

cp = X *-\r — rj|. 

The introduction of the quantity t u \ whicn has no clear physical 
meaning. is necessary only because otherwise we cannot formulate the 
initial electrostatic problem. 

t Since the macroscopic volume average of any quantity in a 
homogeneous system is a constant, then t also cannot depend on the 
value of the coordirate r. The latter is hi^;ily essential, since in the 
literature attempts hav-' repeatedly been undertaken to treat t as a 
function of the distance between ions. This, of course, is false. 



where r„ = r, — r,. Here, in calculating integrals such as 



v 



we have omitted the portion of them that is localized i a 
sphere of radius r 0 centered at r„ since the latter depends 
only on the nature of t he charge distribution within the 
ion. Hence, it is of no interest. 

Averaging (2.2) over the ensemble at constant volume 
V of the system, we obtain 







0de\fJL V ***** 

e d@J \ 2 — ■ er u 






^(r)inr-dr} = (l +? %■ \ ^-<Pr. (2.3) 



This expression represents the electrostatic component of 
the Internal energy of the solution.* It differs from the 
usual formula (1.11) for U, first, by containing the factor 
1 + (0/eXdc/d©), which describes the change in the 
entropy of the folvent in the field of the ions, and second, 
in having the additional term 



* 

characterizing the intrinsic energy of the ions. Both these 
corrections make an appreciable contribution to the 
increment of the internal energy of a solution of finite 
concentration with respect to an infinitely-dilute solution, 
for which t = <•, and 27* = I, 



Avtf <e,) = £/ <e,) (v) - U (tl) (v = 0) 
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e o r o 



(2.4) 



In the special case in which 27* is given by Eq. (1.40), 
A ,£/<*" has the form 



Av/ 



•(H) 



N_ j f . 0 dF_\ j k*r* k-r- 

- I v. F <)Qj I rr 0 F (r„ . r y ) 

_© 1*0 \ *:<■’_ 1 

\ F„ dSj f „r„ f 



1 



(2.5) 



‘The initial formula (2.1) for l/***" can represent either the internal 
energy, or the enthalpy, depending on whether the expression 
( I /4 w)*E dE was integrated at constant volume or at constant pressure, 
in distinction, the formula (2.3) always represents the internal energy, 
since the condition V = const was made explicitly in the averaging. 
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If the dielectric constant r(0, v) = const, then (2.5) 
coincides with the Debye formula (1.12). However, in the 
general case in which t = «(© , v), they differ quite 
substantially. 

II. When the quantity A .U has been found, it is not 
hard to calculate the increment in the free energy of the 
system A/ = F{v) — F ( = 0) by the Gibbs-Helmholtz 
formula U = — ©<<^d©XF/©). Indeed, upon substituting 
into it the A ,U from (2.5) and integrating the expression 
obtained between the limits © = © , and © = © 2 , we find 
that 

M-Y-i 

A«.Y | v(y. 7 .»)~ 1 l»(l x) x * x- | ] , 

(2.6) 

where Ag( } = { fe-e, ~ { le-e,* Th* subscript 0 
applied to x implies that in x = kw/e ©r 0 the quantity c 
= t(v) is taken at v = 0, and is equal to to = e(0). We 
note now that the free energy F of the system is a 
single-valued state function. Hence, the value of A,F(0) 
at 0 = © , cannot depend on the arbitrary temperature 
© = © 2 . The latter condition is compatible with (2.6) only 
if 



A V F V© jy (y. — x | *" ( * x)-‘ y - v | 

/(«>-/(<»)} X/- 4 '" xr"'. (2.7) 

Here A r F (S * = N 0[/(/i) — 7(0)1 and the function An), 
which is an arbitrary integration constant (since A An) = 
0), can depend only on the relative concentration n = 
N/N„ where N, is the total number of solvent molecules 
and N is the total number of particles of the solute. Since 
the expression enclosed within curly brackets in (2.6) is 
valid only at concentrations not exceeding one mole/liter 
(since always <c < 1), then we can assume to a higher 
degree of accuracy that An) — AO) = Bn. 

If we assume in (2.7) that x = X<* then the term Vi(x — 
Xo) first introduced by Huckelm vanishes, and we arrive 
anew at the ordinary Debye expression for the free energy 
increment. We emphasize that this applies particularly to 
the free energy A ,F, but not at all to the thermodynamic 
potential A, = A/* + FA,F, as has often been 
stated.i^i And although the numerical difference 
between Xfr and A f is small (because liquids are 
practically incompressible and A .V ^ 0), it is very 
important, since it means that the chemical potential of 
the solvent has been calculated incorrectly heretofore (see 
below). 

III. Using the ordinary formulas of thermodynamics, 
we can find from (2.7) the rest of the thermodynamic 
functions of the solution: its entropy, enthalpy, etc. 
However, we shall not do this, but shall restrict ourselves 
to treating only the chemical potentials of the solute n 



and the solvent m<- We must bear in mind in calculating 
them the fact that A F<*« = 7(0 , V, N, «, «o), and in turn, 
that € = <(©, V, N, N,), and t 0 = lim# « <(©, K N, N,). 
Since by definition m = (dF/dA0®.v.iv and m, = 
(dF/dAUe. r .w,we can neglect the dependence of A F, t, and 
Co on © and V in the differentiation. Hence, 

<i\y 



•>\y 


o\F tty 


o\r 


or 0 


oX 


tty 







V's 


Ms M , . , V 

id 


* i \y '• 

■' Vs ) 




«\y 


tt\F tty 


o\y 


or,, 


»x s 


~r~7i X~ 




"A's ’ 



( 2 . 8 ) 



(2.9) 



where ^ = j*io» + © In n and m.®’ + © In (1 — n) s m.®’ — 
©.are the chemical potentials of the solute and the solvent 
for the case of an ideal solution, dto/dN = lim* « 
(dt/dN), and dtJdN, = lim v -o (dc/dN,). In taking the 
derivative in (2.9), we must bear in mind the fact that 
only A/**' = ®B(N l /N,) depends directly on N, in (2.7); 
on the other hand, the electrostatic component A.F 111 = 
B(N 2 /N,) depends directly on the electrostatic component 
A /■«*•' of the free energy depends only indirectly on N, via 
( and Co.* Substituting (2.7) into (2.8) and (2.9), we obtain 
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where n = ^ + © In/, and /is the activity coefficient; /i< 
= hj * i - ©„, and g is the osmotic coefficient. 

By using (2.3), the electrostatic component of the 
chemical potentials can be written in the form 



\u 



(2.13) 
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where r = NJV 0 is the density of the solvent. Thus we 
see that Am* 11 is determined (as calculated per ion) by the 
increment in the total electrostatic energy (the integral 
f v (c£ 2 /8t) dV, minus the portion due to the change in 
the dielectric constant of the solvent (the integral f t \ 



•And of course, via V = Nv + N,v„ where v and v, are the partial 
molar volumes of the solute and solvent, respectively. However, we 
need not take into account the dependence of A,P*"on N, via * ^ V~ l \ 
since V is held constant in the differentiation 
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(£V8x)JV(de/dA0 dV), while M, (d ’ is the electrostatic work 
done in changing the density of the solvent. As we should 
expect, the expression for **.'*" exactly coincides with the 
usual definition of the chemical potential of a dielectric in 
an external field.' 2 *' 

IV. If we assume in (2.10) that e = e„ = const, then we 
get the well-known Debye expression for the activity 
coefficient with Onsager’s correction* 

In / -4--j Z V (2.14) 

whereas the expression (2. 1 1) for the osmotic coefficient g 
does not transform under any conditions into the 
formula used in the theory of electrolyte solutions.' 2021241 
In order to understand what the trouble is here, let us 
examine the course of the arguments usually used in 
deriving the expression for g. First, one assumes at the 
outset that (2.7) gives the increment in the thermodynamic 
potential A,.y~, which is wrong, as we have seen. Then one 
finds Afi, by the formula 
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Here one drops the terms containing dt/dN, and deo/d/V, 
for unknown reasons, although, as we know, (t/«X3«/3t) 
^ I -*'°* The error in this derivation is evident. Even if we 
define m. as (3.78 dNfesj,. then by virtue of the identity r 
= F + PV , P = —dF/dV, the term in n, containing 
d /dV nevertheless vanishes, and we again arrive at Eq. 
(2.9). This is natural, since p = (dF/dN) & „ = 
(dP/'dh 

V. Now we shall show that neglecting the con- 
centration-dependence t = t(N) leads to inner contra- 
dictions in the case of electrolyte solutions of intermediate 
concentrations for which k it 0.1. To do this, we recall 
that when P % 0 = const, the chemical potentials n, of the 
solvent and m of the solute satisfy the Gibbs-Duhem 
relation rt*(dfi/bn*) + (I — «*XdM»/dw*) = 0, where n* 
= N/(N + N.) is the so-called mole fraction. Upon 
transforming from n and it, to In / and g, we obtain 1 — g 
+ (I — n*) In (I — n*Xdg/dnn) + ( d In f/dn*) = 0. 
When N « N„ so that n* ^ n, this leads to the 
well-known Bjerrum relation' 21 ': 

\-f! -n~H-g In/I. (2.15) 
To start with, we shall assume that < = e(N.) is 



•On safer suggested that in (2.14 the term Bn coincides with the 
second virial coefficient ?p/3 of a system of hard spheres. Hence B = 
V*rr>N/VKN,/ N) = 4W/6XN,/ V) = V 0 /4V„ where V 0 = »r 0 V 6 
and V, = V/N,. 



independent of the concentration n — N/N,* In this 
case, substituting (2.10) and (2.11) into (2.15) gives 



Such an equality is impossible, since the left-hand side of 
(2.16) is a constant, while the right-hand side 
substantially depends on v by way of k ~~ VT. This 
contradiction can be removed only by assuming that «(v) 

const. In fact, if < = «(t, v), then substitution of (2.10) 
and (2.11) into Bjerrum’ s relation (2.15) transforms the 
latter into a partial differential equation in the unknown 
function <. By solving it, we can find e(r, v), apart from 
an arbitrary function <«» = v). 

Equation (2.16) implies that (r/eXde/dr) = 1/3 for 
extremely dilute solutions, for which «c — * 0. Upon 
substitution into (2. 1 1 ), this gives 1 — g = '/* x*- Since 
this expression, which coincides with the well-known 
Debye-HQckel formula, agrees well with the results of 
measurements,'* 21 - 24 ' we can consider that in this case Eq. 
(2.1 1) admits of experimental verification. However, this 
becomes impossible at higher concentrations, since (2.11) 
contains the unknown derivatives de(r, v)/dr. In order to 
remove this difficulty, we can determine the osmotic 
coefficient g by integrating (2.15): 

*-«— <2i7 > 
0 

where In/ is given by (2.10). The expression thus 
obtained, which contains only e and de/dN , will evidently 
be exact, since it is transformed into an identity upon 
substitution of (2.10) and (2.17) into (2.15). 

12. CoaparisoM with Eiptri f tw 

I. Before going on to comparing directly the formulas 
derived above with the experimental data, we shall make 
some preliminary remarks. 

First, the final formula (1.40) for the binary 
distribution formula satisfies the system of equations 
(1.35), (1.-36) better as x and k becomes smaller. If we 
assume that the system (1.35), (1.36) ensures the 
necessary accuracy, then the discrepancy between the 
thermodynamic and experimental curves must sys- 
tematically increase with increase of these parameters. 
Here it should become very large for * Z 1 and x ~ 3 
(since « s 1 and x = 3 determine the limits of 
applicability of the derived solution). In essence, only this 
sort of agreement can serve to prove the correctness of 
the initial system (1.35), (1.36). However, if the error of 
the theory proves to be independent of the values of x 
and x, this might mean only that it is due to some causes 
other than a systematic error in the derived solution. 



•Since the constant < results from averaging over the N, solvent 
molecules, it cannot fail to depend on N, (or more exactly I, on the 
density r = NJ F 0 ). 




[22 



STATISTICAL THEORY OF ELECTROLYTE SOLUTIONS 



183 



Second, in the Debye model of an electrolyte solution, 
all the forces of non-Coulombic origin are subsumed 
under an intrinsic volume v 0 = xr 0 V 6 occupied by the 
particles of the solute. However, real systems can 
manifest also other forces of a similar sort,* not explicitly 
considered in this model. In part, the effect of these 
additional (with respect to the chosen model) forces can 
always be compensated by varying the value r 0 with 
respect to the crystallographic diameter r c of the ions. 
Hence, in experimentally testing the theoretical formulas, 
r 0 must be treated as a certain adjustable parameter 
effectively characterizing all forces of non-Coulombic 
origin acting in the system. 

Third, any theory taking the solvent in to account only 
i terms of its dielectric constant involves the arbitrary 
function t = t(v, 0). Generally speaking, the form of the 
latter can be determined by independent measurements 
and appropriate calculation. However, such measurements 
are practically lacking at present, rendering impossible an 
all-sided test of the theory. The concentration dependence 
< = «(*>) has been determined only at 25 °C with more or 
less accuracy for aqueous solutions of certain alkali 
halides. Here it turned out that* 25 ’ 

r»(l -f,v). (2.18) 

where v is measured in moles/liter, t 0 = 78.5 is the 
dielectric constant of pure water, and the mean value of e, 
for solutions Qf LiCl, NaCl, KC1, RbCl, NaBr, and Nal 
is 0.20 ± 0.03 (there are no data for the other alkali 
halides). Using (2.18), we can calculate only the activity 
coefficient (2.10)* 



In / 



y*< l v 

(I *|V>- ( I I r,v y*r t , l v) 



7* 

2r 0 ' 



F.V 



(I r.v)2 






(2.19) 



where at 25 C, X o* = 7.134, *,* = 0.3289, r 0 is 
measured in A, and 



X 



/<• 7* . 

1 f,v r„(| F jV) 



X 



x .1 



l I »i v 



o 



l v 

\ 1 'I'' 



( 2 . 20 ) 



The remaining formulas of the previous section contain 
the unknown derivatives 9«/9 r or 9c/ 90, preventing any 
comparison of them with the experimental data. 



*E.g., repulsive forces due to overlap of the hydration shells of the 
ions, or ion-dipole attractive forces (see below). 

•Here, according to Onsager,i»> the product Bn in (2.10) has been 
assumed equal to ftp: the more general case in which Bn ft %p has also 
been discussed in wi. 



II. In order to show the nature of the influence of the 
concentration-dependence of the dielectric constant, we 
shall begin by discussing the very simple case c = c(v) = 
const (i.e., c, = 0). Here, using (2.14) and (2.17), we get 
the ordinary Debye expression for the osmotic coefficient 

g- 

1 g ~2yj 1 1 *-(l X) 1 2 111 (I x)l -z (2.21) 

As was shown above, this gives qualitative results hardly 
differing from Eq. ( 1 .43), which was found firectly from 
the expression ( 1 .40) for the binary distribution function 
by using (1.19). 

Equation (2.21) (or equivalently, Eq. (1.43)) was 
compared with experiment as follows: the theoretical 
curve g = g(r 0 , v) was fitted to the experimental curve at v 
= 0.1 moles/liter by proper choice of r 0 . Then, assuming 
r„ = const, the difference Ag = g ltm - g ap was found at v 
= 0.5 and 1 .0 moles/liter. The values of go, were taken 
from the tables in < M >. The results of the corresponding 
calculations for aqueous solutions of 15 alkali halides are 
given in Table II.* We see from these data that (2.21) 
gives an error that on the average is fourteen times the 
experimental error 8 = 0.001, even at v = 0.5 
moles/liter. At v = 1.0 moles/liter, the discrepancy 
increases to a factor of thirty. Furthermore, there is no 
correlation whatsoever between Ag and the parameter \ 
= c 2 /<0r o . Thus, the experimental test of Eq. (2.21) has 
shown it unsuitable for describing electrolyte solutions of 
intermediate concentrations. However, we should have 
expected this (we recall that in deriving (2.21) the 
assumption was made that t(v) = const, which is 
inadmissible). 

III. It is of definite interest to compare Eq. (2.21) with 
other currently known expressions for the osmotic 
coefficient. As was shown above, the most rigorous of 
these is the Mayer-Haga formula (III). However, the 
series in (III) converges rapidly enough only when v ~ 
0.05 moles/liter, i.e., in the concentration range in which 
the correction to the limiting Debye-Huckel law is so 
small that it can’t be determined experimentally to a 
sufficient degree of accuracy. Hence there is no point in 
comparing (2.21) with the Mayer-Haga formula. 

The next in order of reliability is the Tyablikov- 
Tolmachev formula (1.42) for the distribution function. 
At low concentrations it gives results agreeing with the 
Mayer-Haga formula. Of course, (1.42) is approximate at 
higher concentrations, since it does not include all terms 
of the same order of smallness. However, as was shown 
in <*», it can be derived quite rigorously by selective 
summation of a definite class of diagrams. If we assume 
that the rest of the diagrams contribute only little to the 
final expression, then we can consider the Tyablikov- 



•All the calculations were performed by Yu. M. Kessler, for which 
the author is highly grateful. 
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Table II 
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Name 


V - 


0.1 moles/liter 


V 


0.5 moles/liter 


V = 


-1.0 moles/liter 
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r 0 (A> 




’ ? theor *exp 




A * " ^theor *exp 
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A 


Eq. 


Eq* 


A 


Eq. 


Eq. 


A 


Eq. 


Eq. 


; lyte 




(2.20 


12.211 




C2-23) 


(J21) 




(2-^3) 


(221) 


; Lil 


0.04 


4.43 


4.54 


0.23 


0.(H*9 


0.009 


0.41 


0.035 


0.037 


• l.iltr 


0.04 


4.04 


4. 10 


0.2(» 


0.014 


0.011 


0.45 


0.(»25 


0.019 


1 l.iCI 

i 


0.04 


2.S(i 


.{.‘14 


0,27 


0.007 


((.000 


0.47 


0.010 


0.004 


Nal 


0.04 


3.81 


2.89 


(».27 


0.014 


0.007 


1 

0.48 


1 

0.032 


0.020 


i Nal Sr 


o.o4 


.'(.iso 


:ui5 


0.2 H 


0.018 


0.012 


0.50 ! 


! 0.044 


0.027 


1 Nai:i 


0.o4 | 


| 3.53 


.1 . 57 


0.2H 


0.024 


0.UI8 


0.51 j 

! 


| 0.058 


0.039 


KI 


0.04 


3.53 


.«. 5(i 


0.28 


0.028 


0.022 


0.51 


o,(h;7 


0.048 I 


KlSr 


0.04 


! :s.:s2 




o.:so 


0.027 


0.020 


0.54 


O.IM52 


0.042 


KCl 

1 


o.o4 


j :s.2T 


::.27 


(1.30 


0.028 


0.019 


0.55 


0.1K>7 


0.044 


! ItliCI 


1 

j o o:< 


2.07 


2.07 


0.31 


0.025 


1 

{ 0.010 


0.58 


0.050 


0.033 


• Ittil Sr 


. o.o:S 


2.110 


2. 9S 


0.32 


0.020 


I 0.01! 


0.00 


0.051 


0.028 1 


! Rbl 

1 


j 0.1,5 


2.114 

« 


2.112 


0.32 


0.020 


j o.oio 


0.01 


| 0.048 


0.025 j 
1 


CsCI 


0.02 


2.70 


2.75 


0.34 


0.020 


1 0.010 


1 

0.03 


0.050 


0.027 


(SI Sr 


0.02 


2-74 


2.70 


0.34 


0.027 


I 0.018 


0.04 


o.oSO 


0.032 


Csl 


0.02 


2.70 


2.0*4 


0.34 


0.027 


j 0.017 


0.05 


0.050 


0.031 


i 


r u 


r u 


0.025.4 


! 


r 0.021; 


0.014 


1 


\g 0.048 


; 0.O30 J 



Tolmachev formula to be a sort of extrapolation valid for 
any concentration. We cannot determine directly the size 
of the error that we make in such an extrapolation, since 
this would require that we find an exact expression for 
the binary distribution function. However, we can 
estimate it indirectly to the same degree of accuracy to 
which, e.g., the condition of neutrality (1.5) is obeyed. By 
substituting (1.42) into it and transforming to dimen- 
sionless variables, we obtain 

A-*- * ! S sh [K-T-]'* rf ' < 2 22 > 
1 

Table II gives values of A calculated by this formula for * 
= 0.3, 0.7, and 1 .0 (this corresponds approximately to v 
= 0.1, 0.5, and 1.0 moles/liter).* Since when the 
Tyabilikov-Tolmachev formula is accurate enough: 
however, on going to more concentrated solutions, the 
value of v = 0.1 moles/liter, A is very small,t we can 
assume that at this concentration the Tyablikov-Tolmachev 
formula is accurate enough: however, on going to more 
concentrated solutions, the value of A increases sharply, 
and hence the accuracy of (2.22) must decline. 

Table II gives the values of r 0 and A g calculated by the 
formula t 

bb 

1 - g if l sh [ X -~j— ] ' dt ~ ch f (2 23) 



•These data are taken from ,JI *. 

tWe recall that A must be identically zero in an exact theory. 
$The values of Ar are taken from (1)| . 



obtained upon substituting (1.42) into the expression 
(1.19) for the pressure. The value of r 0 was determined by 
fitting the experimental and theoretical curves at v — 0. 1 
moles/liter. As we see from the data given in the table, 
practically coincide (the arithmetic-mean value of the 
difference Ar 0 is only 0.025 A). However, Eq. (2.21) 
proves to be more accurate at high concentrations than 
the Tyablikov-Tolmachev formula (2.23), since at v = 0.5 
and 1.0 moles/liter the value of Ag calculated by Eq. 
(2.23) is about 1 .5 times as great as the value found by 
Eq. (2.21). 

IV. We shall now proceed to test Eq. (2.21), which was 
derived under the assumption that € = t(v) = var. Table 
III gives the values of r 0 calculated for aqueous solutions 
of 15 alkali halides by fitting the theoretical and 
experimental values of log /at v — 0.1 moles/liter.* The 
effective diameters r d , of the ions are given in the same 
table for comparison, as found from electrical-conductivity 
data.i»i if we consider that in the former case the ionic 
diameters are found by studying the equilibrium 
properties of the solutions, but in the latter case, by their 
kinetic properties, we must acknowledge the agreement 
of the two quantities to be strikingly good. 

A general graph of the relation of log / to the 
concentration v of the electrolyte is given in Fig. 2. We 
see that it shows good qualitative agreement between the 
theoretical and experimental curves; however, the 
quantitative agreement, which is very good for LiCl, 



•The values of /are taken from < M l 
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Table III 



No. 


Name 

of 

salt 


r 0 

(A) 


r eff 

(A) 


Xo 

7.1.14 

r<* 


r .< 

<A> 


\>n 

- r .c 


1 


Ltl 


4.48 




1.59 


2.88 


1.60 


2 


Li Hr 


4.00 


— 


1.78 


2.64 


1.36 


3 


LiCI 


3.89 


3.27 ±0.13 


1.84 


2.49 


1.40 


4 


Nal 


3.69 


3.81 ±0.26 


1.94 


.3.18 


0.51 


5 


NaBr 


3.45 


3.45 


2.07 


2,94 


0.51 


6 


NaCI 


3.30 


3. 28 ±0.22 


2.16 


2.79 


0.51 


7 


KI 


3.23 


3.67±0. 19 


2.21 


3.53 


—0.30 


8 


KHr 


3.00 


3.301-0.08 


2.37 


3.29 


-0,29 


9 


KCI 


2.90 


3.06 r0. 16 


2 . 46 


3.14 


—0.24 


10 


KbCI 


2.59 


2.90 


2.76 


3.30 


-0.71 


11 


Kbltr 


2.54 


— 


2.81 


3.45 


-0.94 


12 


Rbl 


2.45 


— 


2.91 


3.69 | 


1 -1.24 


13 


CsCI 


2.04 


2. HI 


3.49 


3.46 


- 1.42 


14 


Cslir 


2.04 


— 


3.49 


3.61 


— 1.57 


15 


Csl 


2.04 


— 


3.49 


3.85 


1.81 



gradually deteriorates in the series Li > Na > K > Rb 
> Cs. Here we observe a distinct correlation between the 
error of the theory A log / = \ogf lhmt — log and the 
values of the parameters k and x (Fig. 3). The latter is a 
convincing proof of the idea that the disagreement 
between theory and experiment in this case is due to a 
systematic inaccuracy of the obtained solution. However, 
this inaccuracy is not so great, since the relative error 
A log / of the theory nowhere exceeds 13 percent, even at 
v = 1 mole/liter.* 

V. As was noted above, the parameter r 0 effectively 
characterizes the value of all the forces acting in the real 
solution. At the same time, the mean crystallographic 
diameter r e = where the r/p are the 

crystallographic diameters of the anion and cation, 
respectively, defines the value of only those forces 
involving the excluded volume occupied by the ions 
themselves. Hence it is natural to assume that the 
difference A r 0 = r 0 — r c is proportional to the additional 
forces acting in the real solution, and not taken into 
account explicitly in the charged-sphere model. And if 
this is so, then some correlative dependence should exist 
between A r 0 and the parameters determining the strength 
of the additional forces. 

One of the sources of the additional forces in a real 
solution could be the hydration of the ions.i»- M i Two 
particles of the solute can approach closely only after 
their hydration shells have been disrupted, which takes a 
certain expenditure of energy. Then evidently, hydration 
must give rise to additional repulsive forces that can be 
compensated only by an increase in r 0 with respect to r c . 



•The relative error is 13 percent only for CsCI. It drops to aelO 
percent for RbCI, KCI, and NaCI, and is still less for LiCI. 




FIG. 2. Concentration-dependence of the activity coefficients of the 
alkali chlorides. experiment, theory. 



Hence, correlation between the value of A r 0 and the 
hydration energy of ions must be observed most 
clearly for those salts for which A r 0 ^ 0. 

Now we recall that hydration is due to the influence of 
the electric field of the ions on the solvent molecules, 
which have a permanent dipole moment. Other 
conditions being equal, the energy of this interaction is 
proportional to the electric field intensity at the center of 
the water molecules. The latter in turn is inversely 
proportional to the square of the distance /•,_„<*> = [Vi(r f <±> 
+ r Hl o)] between the centers of the ion and of the solvent 




FIG. 3. Relation of the absolute error A log/ = lot f.\ — I o*/„p 

to the parameter x- 
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FIG. 4. Relation of A r 0 = r 0 - r c to the hydration energy of the 
cations {/<*»*> ~ +>) -2 . 

molecule (r H2 o = 2.90 A is the diameter of a water 
molecule). Hence we see that the anions, which have a 
large diameter r f <->, are almost unhydrated, since their 
hydration energy (/<*>*> ~ (r,_j->)- 2 is very small.* 
However, it is considerably larger for the cations, and 
furthermore, for them varies over the rather large 
range (0.10 | ^ 0.20). Hence A r 0 should be 

distinctly correlated only with the hydration energy of the 
cations, as is well confirmed by the factual data. Indeed, 
Fig. 4 shows that the value of Ar 0 increases practically 
linearly with increasing for salts of Li + , Na + , and 
K + . On the other hand, the salts of Rb + and Cs + , for 
which A r 0 is negative and large, fall off the curve. 

We shall now examine what happens to the parameter 
r 0 in the series MCI, MBr, MI, where M is one of the five 
metal cations (i.e., Li, Na, K, Rb, or Cs). Since the 
cations are hydrated but the anions practically not, then 
when the Cl~ anion in the salt MCI is replaced by Br* or 
1“ , the increase in r 0 must be close to the increase in the 
radius (rather than the diameter) of the anion In 
other words, the difference V^[(r 0 ) M i — foWil must be 
approximately equal to V4[r,*s,> — r f < ci s,>] = 0.15 A. We see 
from the data given in Table IV that this requirement is 
satisfied quite well for the Li, Na, and K cations, but not 
for Rb and Cs. All of this leads to the concl ision that the 
additional forces in Rb and Cs salts do not arise from the 
hydration effect, but from some other cause. 

Since the value of A r 0 is negative for Rb and Cs salts, 
we might naturally assume that for them the additional 
forces are attractive. For example, the latter might be due 
to ion-dipole interaction of two solute particles. When 
neither of the particles possesses a permanent dipole 



Table IV 



M 


Li 


Na 

1 


K 

1 


Kl» 


Cs 


y(( r o)MI -( r o)MCll 


0.59 


0.39 


0.33 


- 0.14 


0.00 


'2 l( r o)MBr ( r o)MCI 1 


0.11 

i 


0.15 


0.10 


—0.04 


0.00 



moment, the energy of such an interaction is a£ 2 2 , where 
a, is the polarizability of the first ion, while £ 2 is the 
electric field intensity caused by the second ion at the 
center or the first ion. When the Cl~ ion in a salt MCI is 
replaced by I~ or Br", the value of £ M remains 
practically constant (since E„ ~ l/r f 2 , while r c varies only 
from 3.30 to 3.85 A for the Cs and Rb salts; see Table 
III). However, the polarizability of the Cl~, Br~, and I~ 
ions, being proportional to their volumes (i.e., to (r,<->)’. 
On the other hand, Ar 0 should not depend on (r f <-») J at all 
for the salts of Li, Na, and K, which show practically no 
polarization interaction. This is precisely what is actually 
observed (Fig. 5). 

CONCLUSION 

We shall briefly formulate the results obtained above. 

(A) Depending on their structures, all electrolyte 
solutions can be divided into three groups: dilute 
solutions, solutions of intermediate concentrations, and 
concentrated solutions. The distinguishing feature of the 
first group is the high degree of “collectivity,” such that 
the Debye sphere contains a large number of ions at one 
time. For the second group, the charge of the “central” 
ion is screened by only one counterion, leading to 
formation of neutral quasimolecules in the system. 
Finally, concentrated electrolyte solutions are characterized 




•In the series Cl“ , Br“ , I" , the value of (r ( _ m «->)- 2 ranges only from 
0.078 to 0.098. 



FIG. 5. Relation of Ar 0 = r 0 - r c to the volumes of the anions V - ’ 
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by a structure of the ionic subsystem greatly reminiscent 
of the structure of ordinary liquids in its close-range 
order and coordination spheres. 

(B) Since the absolute values of the density of the ionic 
subsystem in dilute solutions and solutions of intermediate 
concentrations are small enough, one can describe the 
solvent in them in terms of the dielectric constant « = 
e(0 , v) alone. Here one can neglect in dilute solutions the 
dependence of e on the concentration v (but not on the 
temperature 0!). However, in solutions of intermediate 
concentrations, one must always take into account the 
fact that e = e(v). One can’t introduce the dielectric 
constant at all in concentrated solutions, since the 
discrete structure of the solvent plays the important role 
in them. 

(C) One can use methods of expansion in series in a 
small parameter in constructing a theory of dilute 
solutions, but these methods are not suitable for 
describing solutions of intermediate concentrations or 
concentrated systems, since then the corresponding series 
diverge. Hence, the theory of the latter can be based only 
on methods analogous to those used in the theory of 
ordinary liquids. However, since the density of the ionic 
subsystem is small in solutions of intermediate con- 
centrations, one can use the linear equations of the theory 
of liquids to describe them. The non-linear effects must 
be taken into account for concentrated solutions. 

(D) One can find a solution of the linear equations of 
liquid theory in analytic form for solutions of intermediate 
concentrations. Unfortunately, it is applicable only to 
aqueous* solutions of univalent electrolytes for v < 1 
mole/liter. This solution unambiguously determines the 
free energy of the ionic subsystem of the solution, and 
here it turns out that the derived expression for the 
internal energy exactly coincides with that previously 
found by Debye and Huckel. 

(E) Heretofore, in constructing the thermodynamics of 
electrolyte solutions of intermediate concentrations, one 
has had to take account systematically of the fact that e 
= t(v) for them. Introduction of the concentration- 
dependence into the expression for the characteristic 
functions of the system permits one to eliminate not only 
the inner contradictions existing in the thermodynamics 
of electrolyte solutions, but also to gain satisfactory 
agreement of theory with experiment. 

(F) An analysis of the values of the single adjustable 
parameter of the theory (the ionic diameter r 0 ) permits 
one to establish the fact that in some solutions repulsive 
forces act between the ions, due to overlap of their 
hydration shells (hydration-type systems), while attractive 
forces act in others, due to the mutual polarizability of 
the ions (polarization-type systems).* 

•One usually speaks of positive and negative hydration. However, it 
seems to us that the term “hydration- and polarization-type 
electrolytes" corresponds better to the physical nature of the 
phenomenon. 
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